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Abstract - As part of the NASA’s DESDynl mission, global 
elevation profiles from contiguous 25 m footprint Lidar 
measurements will be made. Here we present results of a 
performance simulation of a single pass of the multi-beam 
Lidar instrument over uplifted marine terraces in southern 
Alaska. The significance of the Lidar simulations is that 
surface topography would be captured at sufficient resolution 
for mapping uplifted terraces features but it will be hard to 
discern l-2m topographic change over features less than tens of 
meters in width. Since Lidar would penetrate most vegetation, 
the accurate bald Earth elevation profiles will give new 
elevation information beyond the standard 30-m DEM. 

I. Introduction 

Although the primary objective of the DESDynl-Lidar is 
the recovery of vegetation structure for estimation of above 
ground biomass and for ecosystem health, the 25 m 
footprint, contiguous multi-beam data [1] will provide 
important new information on surface topography and 
topographic change. Over the last decade we have 
conducted a number of pre- and post-launch studies that 
support the rationale for the Lidar simulations reported here. 
For example, prior to the launch of the Ice, Cloud and land 
Elevation laser altimeter satellite (ICESat, single laser beam 
with 55-80 m footprints), we used NASA’s Airborne 
Topographic Mapper [2] (ATM) to estimate how effective 
individual ICESat footprint data from exact repeat satellite 
tracks would be in recovering ice elevation change over the 
rapidly retreating Jakobshavn Glacier. These simulations 
enabled us to identify an approach to successfully use 
ICESat plus InSAR derived digital elevation models (DEM) 
for ice elevation change on mountain glaciers [3,4] Also, 
we have used repeat Laser Vegetation Imaging Sensor 
(LVIS) measurements with full-waveform, 20-m contiguous 
footprints, to quantify pyroclastic and lava flow of a 
stratovolcano in Costa Rica[5]. In another study, near 
Barrow, Alaska, we introduced high quality ICESat derived 
positions into processing of ERS-1/2 data to improve the 
vertical resolution of the InSAR derived DEM by a factor of 
three [6], An extension of this approach to fuse the 
moderate resolution multi-beam DESDynl-Lidar into 
processing of SAR data for derivation of DEM for a wider 
range of topographic slope and roughness could be 
especially useful outside SRTM coverage. 

The complex plate boundary in southern coastal Alaska 
provides an excellent setting for testing DESDynl 
capabilities to recover fundamental parameters of 


seismotectonic processes. Also, aircraft acquired Lidar and 
SAR data have been acquired in this region in the last 
decade that can be used for simulations. Here we use the 
National Center for Airborne Laser Mapping lidar 
(NCALM) data acquired over marine terraces west of Icy 
Bay (“Sullivan data”) in southern Alaska in September 2005 
as part of the NSF STEEP study (www.ig.utexas.edu/steep ). 
In this paper, we report on the simulation studies of the 
multi-beam Lidar instrument performance over uplifted 
marine terraces, active faults and folds, and possible uplift 
associated with the 1899 Yakataga seismic event (M=8 ) in 
southern, coastal Alaska. 

II. Approach 

In this performance study we created synthetic waveforms 
by extracting the NCALM point cloud data along one set of 
the DESDynl-Lidar calculated ground tracks (Fig. 1). We 
used a single pass of DESDynl-Lidar satellite ground tracks 
with 5 laser beams and across track spacing of -850 m to 
geolocate the simulated waveforms. Each ground track had 
individual laser returns spaced 30 m along track (Fig. 1). 



Fig. 1. NCALM-derived topographic map of study region showing the 
marine terrace surfaces (I-iV) and the DESDynl-Lidar (NE-SW) ground 
tracks. The offset between terrace surface II and III is frequently masked 
by alluvial deposits. A subtle offset of 1-2 m between terrace IV and III 
may be due to coseismic displacement associated with the 1899 M=8 
earthquake. The northern portion of profile #1 crosses the Little River and 
is not included. Inset: outline of Alaska with study area. 
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By accounting for sensor characteristics such as footprint 
size and the temporal and spatial distribution of laser 
energy, we generated pseudo-DESDynl waveforms from the 
last return point cloud data (Fig. 2) using the approach 
developed earlier for aircraft LVIS data [7], Elevation 
measurements are extracted using information obtained in 
the waveform only. The lowest and highest elevations 
within each footprint correspond to the lowest and highest 
detected signals above a predefined threshold. The centroid 
of the waveform has also been extracted to investigate the 
precision and accuracy of this as an elevation measurement 
over the terrace surfaces. These waveform derived metrics 
have been compared to corresponding small-footprint 
elevations to assess their precision and accuracy. 


Sullivan fault is marked by a prominent escarpment on the 
mountain above the highest terrace level; this feature is 
more evident to the west of the region shown in Fig. 1. The 
thick coastal forest and brush mask important aspects of the 
geomorphology on aerial photographs, standard DEM, and 
C-band synthetic aperture radar (SAR) images, a problem 
that is largely overcome by the NCALM Lidar data [9], An 
advantage of even the medium footprint DESDynl-Lidar 
data over the standard existing 30 m DEM is the ability to 
see through the vegetation to discern the individual terrace 
heights, features of the escarpment associated with the 
Sullivan fault, and to map out the morphology of the 
Sullivan anticline. Additionally, the DESDynl-Lidar 
derived elevation profiles will be provided in a common 
global reference frame. This will facilitate fusion of this 
data with other DEM and satellite images. 



Fig. 2. WGS84 elevation versus waveform amplitude for a simulated 
waveform over a footprint located along Track 1 near the transition from 
near-coastal terrace level IV to the next terrace level. The waveform has 
been simulated assuming DESDynl-like parameters for the observations 
(i.e., laser pulse width (9ns), laser footprint diameter (25m), waveform 
vertical resolution (2ns)). Random noise (2 counts) was also included to 
simulate background noise present in the measurement. The range between 
the lowest and highest elevation could be due to a variable ground surface 
or NCALM Lidar returns from low-lying brush (higher curve). 


III. Geodetic imaging of seismotectonic structures 

IN SOUTHERN ALASKA 

Between Icy Bay and Cape Yakataga there are four 
marine terraces that range in elevation from approximately 1 
m to 50 m located on the southern flank of the Sullivan 
anticline, and in the footwall of the Sullivan thrust fault 8 . 
These features trend roughly east-west, a nearly optimum 
orientation for future roughly N-S DESDynl-Lidar ground 
tracks (Fig. 1). The terraces contain evidence for rapid 
uplift during large to great magnitude earthquakes [8], 
Terrace surfaces are hundreds of meters to over a kilometer 
in width, partly buried by alluvial fans, and separated by 
remnant sea cliffs that mark the former shorelines. The 


A. Centroid elevations along five DESDynl-Lidar ground 
tracks 

The centroid elevation for each footprint has been used to 
create profiles along the 5 tracks shown in Fig. 1. For 
example, the Track 2 and 3 profiles (from southwest to 
northeast) in Fig. 3 show offset (in red) of the distinct sloped 
surfaces (blue lines) associated with the 4 different terraces. 
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Fig. 3 Elevation profiles along Track 2 (top) and Track 3 (bottom) as a 
function of latitude with estimated offset (red) and terrace surface slope 
(blue). 
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We examined the centroid elevation plotted in Fig. 3 (blue 
diamond) compared to the highest and lowest elevation 
(above the background noise threshold) from the NCALM 
data. The elevation differences over the flat terrace surfaces 
were generally 2-4 m but over rough areas this difference 
could be larger (up to 9+ m). For the fairly flat terrace 
surfaces, the centroid elevation was generally in the middle 
of the upper and lower elevation bounds from the NCALM 
data. As expect for a 25 m footprint, the upper and lower 
bounds transition regions over faults or the terraces were 
large and variable relative to the centroid. 

The large offsets between individual terrace levels as well 
as the slope trends across terraces in a common reference 
frame enable more robust comparisons. For example, the 
terrace offsets are largest in Track 2 and then the offsets 
become more muted to the west, which cross alluvial fans 
that have formed on the lower terrace level following uplift, 
and partly bury the original sea cliffs (Fig. 1 and 3). In 
general, the individual terraces occur at lower elevations on 
the more westerly tracks, this partly caused by a lateral 
decrease in elevation away from the high points, or apices of 
a temporal sequence of alluvial fans that built out from 
Little River onto each terrace level. Only after accounting 
for lateral variations in alluvial fan aggradation on the 
terrace surfaces can we search for evidence of deformation 
and uplift related to earthquake slip or post-glacial uplift 
episodes to the east near Icy Bay. 

The significance of the DESDynl-Lidar simulations is that 
surface topography is captured at sufficient resolution for 
mapping uplifted marine ten-aces, including the alluvial fans 
that built out and buried the marine deposits following 
uplift. Relic sea cliffs at the back of the terraces are also 
resolved as abrupt elevation steps with topographic relief 
dictated by the amount of post uplift burial by alluvial 
processes, and to a lesser degree erosional dissection. 
Although the moderate-scale topographic features are 
captured, uplift on features such as the 1899 berm that are 
just several tens of meters across in width are hard to 
discern in the DESDynl-Lidar elevation profiles. 

B. Future steps in waveform simulation over the terraces 

To simulate the DESDynl-Lidar waveforms from aircraft 

point cloud Lidar data, returns from all surfaces should be 
included as the full waveform representation is an important 
feature for identifying vegetation structure. The NCALM 
data retains the first and last Lidar returns. In the next step 
of our simulation we will include both the first and last 
returns. We will estimate the ground elevation beneath any 
canopy using an automatic search algorithm. As mentioned 
in the Fig. 2 caption, the first (highest) curve may actually 
be reflections from low-lying brush. An advantage of full 
waveform Lidar systems such as LVIS and DESDynl-Lidar 
is that is it easier to associate the various returns with known 
vegetation structure. 

C. Ongoing simulations over fold and fault structures 
covered by glaciers 


In the late summer of 2005, ATM laser swath (-480 m 
width) data [2] was acquired in coastal, southeastern Alaska 
and it included data from the Malaspina and Seward glaciers 
(https://www.atm. wallops.gov) . We used a version of this 
data (ICESSN standard nadir product) with a larger footprint 
size (80 m) to examine the ice characteristics along the 
ATM profile, In the lower Seward the glacier surface 
elevation, along track slope, and within footprint roughness 
and slope were used to identify thrust faults beneath the ice. 
These results are being used to interpret ice velocities 
derived from ALOS PALSAR L-Band InSAR data [9] 

IV. COSEISMIC SLIP RECOVERY SCENARIO FOR A SOUTHERN 
ALASKA EARTHQUAKE 

As summarized in the DESDynl Science Definition 
Document, “in areas where the surface is disrupted over the 
course of a large geophysical event such as coseismic 
displacement and shaking of the surface, the InSAR-derived 
surface deformations measurements cannot be made. It is 
helpful to take advantage of the Lidar to repeat elevation 
transects over the deformed surface to fill in some area of 
disruption.” [10] As summarized above, we simulated 
waveforms along DESDynl-Lidar ground tracks over the 
uplifted region associated earlier great earthquakes with just 
one DESDynl-Lidar pass assuming a 5 laser beam 
configuration. To estimate coseismic displacements 
following a large earthquake, it would useful to use an 
existing pre-earthquake DEM along with the DESDynl 
Lidar over the region. Alternatively, if an off-nadir pointing 
capability exits, repeat targeting of the near-field coseismic 
deformation region would increase the number of transects 
over a short-coseismic interval. This ongoing study is 
enabling us to assess the minimum number of transects 
needed to provide timely, and spatially dense enough, data 
for estimation of earthquake parameters to calculate 
associated stress redistribution and triggered slip on nearby 
faults, inform detailed field studies, and to enable disaster 
managers to assess damage as discussed in the “Use Case 
Scenario” for a Seismic Event in California outlined in the 
Report of the DESDynl Applications Workshop [11], 

V. Summary 

Since DESDynl-Lidar would penetrate most vegetation, 
the accurate bald Earth elevation profiles will give new 
information beyond the standard 30-m DEM. In our 
southern Alaska study region, the major crustal 
seismotectonic features are roughly perpendicular to the 
satellite ground tracks and the accurate elevation profiles are 
enabling process studies of folds, faults, and marine 
terraces. For a major earthquake with on-land surface 
disruption, such as a thrust-slip event in southern Alaska, 
DESDynl-Lidar could provide new and timely near-field 
vertical static displacements when used in combination with 
a pre-earthquake acquired DEM or with some off-nadir 
pointing to previous DESDynl-Lidar transects. We have 
found that these new simulations inform prelaunch 



strategies for maximizing our ability to expand the range of 
solid Earth applications for this medium resolution Lidar as 
part of the DESDynl mission. 
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